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Abstract 

We analytically discuss a multiplicative noise generalization of the 
Kuramoto-Sakaguchi dynamics for an assembly of globally coupled 
phase oscillators. In the mean field limit, the resulting class of in- 
variant measures coincides with a generalized, two parameter family 
of angular von Mises probability distributions which is governed by the 
exit law from the unit disc of a hyperbolic drifted Brownian motion. 
Our dynamics offers a simple yet analytically tractable generalization 
of Kuramoto-Sakaguchi dynamics with two control parameters. We 
derive an exact and very compact relation between the two control 
parameters at the onset of phase oscillators synchronization. 
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1 Introduction and model 

The modeling efforts aiming to analytically explain the phenomenon 
of collective synchronization of phase oscillators are, to a large extent, 
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based on the Kuramoto model [HE] and its mathematical tractable ex- 
tensions [3]. The success of this "toy model" is thoroughly documented 
in the recent review by Acebron et al. [4]. 

Despite a very large body of available literature on coupled phase os- 
cillators, we feel that intimate connections with directional statistics, 
typically describing distributions on the circle, yet remained relatively 
neglected. The main purpose of this note is precisely to highlight this 
intimate connection and explore the fruitful consequences that can be 
derived from it. As a matter of fact, the mean field stationary solu- 
tion to the noisy Kuramoto model proposed by Sakaguchi, (referred 
from now on as the KS model) does on the circle and despite to an 
apparent non-linearity in the drift term, effectively play the role the 
Ornstein-Uhlenbeck process plays on the line (see discussion in Mardia 
[5]). This remark enables to better comprehend why "this little won- 
der", as the KS model is called in [3] and [6], is analytically tractable. 
This also suggests that a "non-flat" extension of the mentioned circular 
Ornstein-Uhenbeck process could be constructed and used to guide us 
to what we shall call here an hyperbolic extension of the KS model. 
Our starting point is a generalized version of the Kuramoto model 
introduced by Sakaguchi [3]. The model consists of N non- linearly 
coupled phase oscillators having identical natural frequencies uj and 
which are perturbed by multiplicative Gaussian White Noise. The 
phase of the i-th oscillator, denoted by 6i, evolves in time according to 




(1) 

The first summation over the oscillation population in eq.JT]) is a de- 
terministic, non-linear coupling with strength K > 0. The second 
summation in eq.JT]) couples the oscillators to random fluctuations 
through the parameter < C < 1. The resulting state- dependent 
noise strength modulates the amplitude of a Gaussian White Noise 
(WGN) source £j's with constant noise strength v2T i.e., 

(fc(t)) = and <&(*)&(*')> = 2T8 t3 5(t - t'). (2) 

The additional coupling parameter C controls the multiplicative char- 
acter of the noise source. Indeed, for C = the additive WGN case is 
reproduced and eq.(H|) therefore coincides with the original Kuramoto- 
Sakaguchi model. 

For C > 0, we are in presence of multiplicative WGN and the under- 
lying stochastic integral will be interpreted here in the Ito sense. 
Observe that the multiplicative noise in eq.JT]) itself depends on inter- 
actions between the oscillators of the assembly. Several multiplicative 
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noise sources for KS dynamics have been previously introduced. In [7], 
the authors discuss a multiplicative noise KS in which the WGN mod- 
ulation affecting each oscillator is determined by its dynamical state 
only, (local noise modulation). A global noise modulation exhibiting 
similarities with our present model is discussed in [8] where the noise 
source perturbs the coupling control parameter itself. The [5] dynam- 
ics differs from our present model by the fact that in the synchronized 
states, (i.e. when 9{ = Vi, k) the dynamics become deterministic - 
a specific feature absent in our model. 

In general multiplicative noise sources are known to give rise to new 
emerging behaviors, the core mechanism for noise induced phase tran- 
sitions [9]. In the sequel, we shall indeed observe that the onset of 
synchronized behaviour will be strongly modified when C > 0. 



As it is usual, we now analyze eqs. ([I]) and (|2j) in the mean field 
approximation reached in the thermodynamic limit A" — > oo. By going 
into a rotating frame, we can get rid of the common eigenfrequency 
and hence, we take from now on u = 0. The collection of oscillators 
can then be characterized by a circular probability density p(6, t) where 
p(9+d9, t)—p(6, t) stands for the fraction of oscillators on the unit circle 
which, at time t, are located at the angular position 9. The circular 
probability density satisfies the non-linear Fokker-Planck equation 



-£ = — — [RK sin(yj - 9)p] + IW(1 + RC cos(</> - 9))p], (3) 



with periodic boundary conditions p{9 + 2tt, t) = p(9, t) and normal- 
ization Jq V p(9,t)d9 = 1. 

In eq.([3|), the quantity R = R(t) G [0, 1] yields a measure for the phase 
coherence of the oscillators and ip = ip(t) represents the average phase 
of the assembly. Both R(t) and tp(t) are used to construct the complex 
order parameter 



Following the original procedure used for the KS model in [T], we 
use R to monitor the cooperative state of the oscillators assembly. 
When R = 1 fully synchronized motion is obtained. The state R = 
characterizes the fully incoherent behavior. Intermediate states where 
< R < 1 indicates that the assembly is in a partially synchronized 
state. 



2 Self-consistent steady solutions 
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Under steady state conditions, a self-consistency condition fixes com- 
patible values for R and tp. This ultimately enables to extract the 
critical coupling strength K c leading to the onset of partially synchro- 
nized dynamical states. 

For the well known additive noise case C = the stationary solution 
to eq. ([3]) is given by the famous von Mises angular distribution 

„k cos(i/>— 9) 

where Iq is the modified Bessel function of the first kind and where we 
introduced the dimensionless factor 

k := RK/T . (6) 

Solving self-consistently eq.(jU together with ([5]) we find the relation: 

R = h(k)/I (k), (7) 

with h(z) the modified Bessel function of the first kind of order 1. 
For the onset of synchronization, i.e. for small but finite R, we may 
expand Eq.([7|) to first order in R. This yields the well known critical 
coupling K c = 2T from where partially synchronized solutions branch 
off the fully incoherent solution represented by the uniform angular 
distribution p(9) = l/2ir. 

Consider now the generalized, multiplicative noise KS dynamics ob- 
tained when C > in eq.([T]). The stationary solution of eq. (J3]) is 
given by the so called hyperbolic von Mises distribution (see [10], [H] 
and the discussion in the next section): 

p{6) = Z" 1 \l + tanh(? ? ) cos(> - 0)1 ° (8) 

with 

K 

tanh(ry) := RC and a := 1 — — — • 

1 O 

The normalization Z in eq.jH]) involves the Legendre function of 
order and degree — a: 

Z = 27rPl°2(cosh(?7)) cosh(r ? ) a . 

In full analogy with the case C = 0, we solve self-consistently eq.([l| 
together with (jSJ) and find the relation (see appendix): 

1 P^coshfo)) 
1 ~ a Pi 2(cosh(r/)) 
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where is the Legendre function of order 1 and degree —a. To 
localize the onset of synchronization, we again expand Eq.© to first 
order in R, (see appendix) and obtain now the C-dependent relation 

K C = T(2 + C). (10) 

The remarkably simple form of eq. (|10p . which is our central result, 
relates the coupling K and the multiplicative noise C constants with 
the WGN variance T. The presence of the extra control parameter C in 
eq. tfTO]) modifies the phase transition diagram for the synchronization 
regimes and offers therefore a natural and very simple generalization 
of the original Kuramoto-Sakaguchi dynamics. 



3 The hyperbolic geometry interpreta- 
tion of the generalized Kuramoto model 

As mentioned before, the KS model stationary probability measure 
coincides with the von Mises distribution and is a cornerstone distri- 
bution in the theory of directional statistics comparable in many ways 
with the normal distribution on the line (see Mardia |12j). Among 
the numerous possibilities to obtain the von Mises distribution eq.jS]), 
let us emphasize here that it describes the exit law from the disc 
x 2 + y 2 < 1 of a planar Brownian motion starting at the origin subject 
to a constant drift vector u of length k. 

Recently, a very elegant contribution of J.-C. Gruet shows how a two- 
parameter hyperbolic von Mises distribution A a (r, du) can be con- 
structed from the exit law of the hyperbolic Brownian motion with 
drift a £ E from the hyperbolic disc of radius r, centered at i, |10j . 
This hyperbolic extension of the von Mises distribution can be written 
as 

(1 + tanh(r) cos(<9)r Q r 
2^P°j(cosh(r)) 

which is obviously equivalent to the centered (i.e., tp = 0) equation 
([8]). This will establish a connection between the generalized KS model 
eq.JT]) and the hyperbolic von Mises distribution. 



Observe that independently of Gruets' work the distribution eq. ffTTj) 
has been recently introduced by Jones and Pewsey as a new family of 
symmetric distributions on the circle including prominent special cases 
such as the von Mises distribution, the wrapped Cauchy distribution 
and others [11]. This family, call them h"(r,9), may be described 
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via the A Q 's by coupling r with a through the curvature — c of the 
hyperbolic space according to the relations r — ► rc and a — ► a/c: 

/£(r, 0) := A a/c (rc, 0), a£l and c E R + . (12) 

Here r > stands for a concentration parameter. The index — c can be 
interpreted as the curvature of the hyperbolic space. As the curvature 
increases to zero, the usual von Mises distribution is recovered, (see 
e-g-, [H] where they consider the special case a = 1): 

ar cos(0) 

lim h%(r,0) = — — — . (13) 
-c/o cK ' 2irl (ar) v ' 

Now it becomes clear, how to interpret the constant C > in eq.flTJ). 
We have i] = RC + 0((RC) 3 ) and a = 1 - K/TC. Hence, small 
values of C > produce the following representation for the stationary 
distribution eq. (j8j) : 

p(9) = h I < /T {R,9). (14) 

We then conclude that — C may indeed be interpreted as the curvature 
of the non-flat manifold underlying the generalized KS model given in 
eq. (HJ). Note that the hyperbolic model converges for C — > together 
with the stationary distribution ([8]) and the critical coupling (fTOl) to 
its flat analogue originally discussed in [3]. 



4 Conclusions 

It is immediate to realize that the stationary measure of the Kuramoto- 
Sakaguchi (KS) dynamics coincides with the famous von Mises direc- 
tional statistics. This distribution on the circle can itself be inter- 
preted as the exit law from a disc of a drifted Brownian motion. An 
hyperbolic extension, i.e. the exit law of an hyperbolic Brownian mo- 
tion from a Poincare disk produces a generalized, two-parameter von 
Mises probability distribution which has very recently been considered 
both from a stochastic and a purely statistical point of views. The 
additional control parameter is intimately connected to the curvature 
of the underlying non-flat manifold. By analyzing the implication of 
these features for the KS dynamics, we have been able to construct a 
simple and natural extension of the original KS model for which, to 
paraphrase [4] and [6], the "little wonder" remains possible. Our gen- 
eralization includes, besides the coupling strength between the phase 
oscillators, a second control parameter, (i.e. the underlying curva- 
ture of the curved manifold), which modulates a multiplicative noise 
source. Ultimately and in full analogy with the original KS model, 
it is straightforward to derive an exact generalized synchronization 
diagram which now depends on two control parameters. 
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5 Appendix 



An efficient way to derive eq.JT} is to use the following series repre- 
sentation for the hyperbolic von Mises distribution valid for a < 1 

m- 

1 r(l-q) P ( Z\cosHr)) 
0(0) = > —. r tttt coslmiw — 6)). (lo) 

Replacing the above series in eq.(@|) and using a trigonometric identity 
together with an orthogonality argument for the cosine functions, we 
find trough elementary integration of cos(#) and cos(#) 2 over [— n, it] 
the equation 

Re** = 1 p -l(<**Hri)) (16) 
1 ~ a P^icosHrj)) 

which is eq.©. 

Eq. (fT0]) follows from the above by carefully expanding the Legendre 
functions and P^ to first order in R. Recall that tanh(r/) = 
RC and hence 77 = \ In ( ) • Using in the following the tables in 
Abramowitz and Stegun [13] we have with 8.11.1: 



J Pl°2(cosh(r ? )) = e^-^F^I l _ a; l; l _ e " 2 ^ (17) 

For P_ a more care is needed because of a singularity in the Gauss 
series for the hypergeometric functions F. Using 8.11.1 we have 

^(eoshW) = Km f(T L_ y F( - 1, -«; 1 - 2„; 1 - (U 

The limit is taken using the fact that (see reflection and duplication 
formula for T-function 6.1.17 and 6.1.18): 

, r(l-2u) , 2 2 ( 1 ^)~ 1 /2 i i 

lim-^ ^ = lim^= T(l-u+-) = — . 19 
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We then have with 15.1.2 in 1131: 



P?>(co S h(,)) = 4e(«-^ a( .-l) ^ _ „-„) V(§,2 _ a;3; , _ 

Putting eas. pTjl and flU} into eq.© and noting that 1 - e~ 2ri = 
we find in first order in R 

a 2RC (1-&) 2RC 

Al + RC 4 1 + tfC 1 ' 

which reduces for small but finite R to K = T(2 + C). 
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